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We describe the implications of cosmic microwave background (CMB) observations and galaxy and cluster 
surveys of large scale structure (LSS) for theories of cosmic structure formation, especially emphasizing the recent 
Boomerang and Maxima CMB balloon experiments. The inflation-based cosmic structure formation paradigm we 
have been operating with for two decades has never been in better shape. Here we primarily focus on a simplified 
inflation parameter set, {w{,,a> ct j m , Qtot, Qa, n s ,rc, erg}. Combining all of the current CMB+LSS data points to the 
remarkable conclusion that the local Hubble patch we can access has little mean curvature (Qtot = 1-08 ± 0.06) 
and the initial fluctuations were nearly scale invariant (n s = 1.03 ± 0.08), both predictions of (non-baroque) 
inflation theory. The baryon density is found to be slightly larger than that preferred by independent Big Bang 
Nucleosynthesis estimates (u b = fl b h 2 = 0.030 ± 0.005 cf. 0.019 ± 0.002). The CDM density is in the expected 
range (oj ct j m = 0.17±0.02). Even stranger is the CMB+LSS evidence that the density of the universe is dominated 
by unclustered energy akin to the cosmological constant (Qa = 0.66 ± 0.06), at the same level as that inferred 
from high redshift supernova observations. We also sketch the CMB+LSS implications for massive neutrinos. 



1. CMB Anisotropies & Distortions 

1.1. What Was Revealed in 1999/2000 

As we penetrate the CMB, we go from a nearly 
perfect blackbody of 2.725 ± 0.002 K ||, to a 
3.372 ± 0.007 mK dipole pattern associated with 
the 300 km s _1 flow of the earth in the CMB, to 
the now familiar motley pattern of anisotropies 
associated with 2 < I < 20 multipoles at the 
30/iK level revealed by COBE at 7° resolution, 
then at higher £ in some 19 other experiments 
(though most with many fewer resolution ele- 
ments than the 600 or so for COBE). The pic- 
ture was dramatically improved this year, as re- 
sults were announced first in summer 99 from 
the ground-based TOCO experiment in Chile 
then in November 1999 from the North Ameri- 
can balloon test flight of Boomerang || , then in 
April 2000 from the first CMB long duration bal- 
loon (LDB) flight, of Boomerang Q, followed in 
May 2000 from the night flight of Maxima [g. 
Boomerang's best resolution was 10', about 40 
times better than that of COBE, with tens of 



thousands of resolution elements. Maxima had 
a similar resolution but covered an order of mag- 
nitude less sky. Both experiments were designed 
to reveal the primary anisotropies of the CMB, 
those which can be calculated using linear per- 
turbation theory. What we see in Fig. |l| are two 
images of soundwave patterns that existed about 
300,000 years after the Big Bang, when the pho- 
tons were freed from the plasma. The visually ev- 
ident structure on degree scales is even more ap- 
parent in the power spectra of the Fourier trans- 
form of the maps (Fig. |^), which show a dom- 
inant (first acoustic) peak, a less prominent (or 
non-existent) second one, and the hint of a third 
one from Maxima. In the following, we describe 
the implications of these results for cosmology. 
Space constraints preclude adequate referencing 
here, but these are given in 

We are only at the beginning of the high pre- 
cision CMB era. HEMT-based interferometers 
are already in place taking data: the VSA (Very 
Small Array) in Tenerife, the CBI (Cosmic Back- 
ground Imager) in Chile, DASI (Degree Angular 
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Figure 1. The Boomerang 150 A GHz bolometer ma] 
(one out of 16) is shown. Of the entire 1800 squar< 
degrees covered (not all of which is shown), only th< 
interior 440 sq. degs. within the central rectangl 
were used in our analysis, i.e., in all less than 59 
of the data. The 124 square degree Maxima-I map 
drawn to scale, is also shown. Bandpowers derive! 
from these maps and their noise correlation matrice 
(Figure Q) were used to derive the parameters givei 
in Table (l|. 
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Scale Interferometer) at the South Pole, where 
the bolometer-based single dish ACBAR exper- 
iment will operate this year. Other LDBs will 
be flying within the next few years: Arkeops, 
Tophat, Beast/Boost; and in 2001, Boomerang 
will fly again, this time concentrating on po- 
larization. As well, MAXIMA will fly as the 
polarization-targeting MAXIPOL. In April 2001, 
NASA will launch the all-sky HEMT-based MAP 
satellite, with 12' resolution. Further downstream 
in 2007, ESA will launch the bolometer+HEMT- 
based Planck satellite, with 5' resolution. 

1.2. Primary &; Secondary Processes 

As we peer back in time, we reach a fuzzy wall 
at redshift z ~ 1100, the photon decoupling "sur- 
face" , through which the Universe passed from 
optically thick to thin to Thomson scattering over 
a physical distance ~ 10 kpc, i.e., ~ 10 Mpc co- 
moving. Prior to this, acoustic wave patterns in 
the tightly-coupled photon-baryon fluid on scales 
below the "sound crossing distance" at decou- 
pling, < 100 kpc physical, < 100 Mpc comoving, 



Figure 2. The Ct (defined in terms of CMB tem- 
perature anisotropy multipoles, (AT /T)i m , as in- 
dicated) grouped in bandpowers for prior-CMB ex- 
periments, including TOCO and the North Ameri- 
can Boomerang test flight, (squares) are contrasted 
with the Boomerang-98 LDB (crosses) and Maxima- 
I (triangles) results in the upper panel. The lower 
panel shows the optimally combined power spectra 
of Boomerang+Maxima+DMR (squares) and con- 
trasts it with that for Boomerang+Maxima+prior- 
CMB (circles), showing that including the prior ex- 
periments does not make a large difference to the re- 
sults. Best-fit models for arbitrary Q to t and for the 
restriction to the flat case are shown in both panels. 



3 



were viscously damped, strongly so on scales be- 
low the ~ 10 kpc thickness over which decoupling 
occurred. After, photons freely-streamed along 
geodesies to us, mapping (through the angular 
diameter distance relation) the post-decoupling 
spatial structures in the temperature to the an- 
gular patterns we observe now as the primary 
CMB anisotropies. The maps are images pro- 
jected through the fuzzy decoupling surface of 
the acoustic waves (photon bunching), the elec- 
tron flow (Doppler effect) and the gravitational 
potential peaks and troughs (" naive" Sachs- Wolfe 
effect) back then. Free-streaming along our (lin- 
early perturbed) past light cone leaves the pat- 
tern largely unaffected, except that temporal evo- 
lution in the gravitational potential wells as the 
photons propagate through them leaves a further 
AT imprint, called the integrated Sachs- Wolfe ef- 
fect. Intense theoretical work over three decades 
has put accurate calculations of this linear cosmo- 
logical radiative transfer on a firm footing, and 
there is a speedy, publicly available and widely 
used code for evaluation of anisotropies in a wide 
variety of cosmological scenarios, "CMBfast" §. 

Secondary anisotropies associated with nonlin- 
ear effects also leave an imprint. The photons 
were weakly-lensed by intervening mass, were 
Thompson-scattered by the flowing gas once it 
became "reionized" at z ~ 20, and were occasion- 
ally Compton-upscattered, with mean fractional 
energy shifts 4T gas /m e c 2 , by nonlinear hot gas 
(thermal Sunyaev-Zeldovich or Compton cooling 
effect); the flowing cluster induces a further ki- 
netic SZ effect. Of these secondary processes, all 
have the same spectral signature of a perturbed 
blackbody, with a A-independent thermodynamic 
temperature perturbation, AT(A, q) — except for 
the thermal SZ effect. 

The SZ spectral distortion is related to the 
Compton y parameter, which is proportional to 
the line of sight integral of the electron pressure: 
AT is negative, — 2y, on the Rayleigh- Jeans side 
and positive on the Wien side, (E 1 /T cm b)y, pass- 
ing through zero at 1378/^m or 218 GHz. The 
COBE FIRAS experiment limits the energy in- 
jection that leads to the thermal SZ pattern to 
be SEcompt on cool 

/E cmb = 4y< 6.0 x 10~ 5 (95% 
CL) [Q. It is routine now to see the (concen- 



trated) fluctuations in y from clusters, though 
the sky-average y is predicted to be < 10 -5 . SZ 
anisotropies have been probed by single dishes, 
the OVRO and BIMA mm arrays, and the Ryle 
interferometer. A number of HEMT-based in- 
terferometers being built are more ambitious: 
AMI (Britain), the JCA (Chicago), AMIBA (Tai- 
wan), MINT (Princeton). As well, other kinds 
of bolometer-based experiments will be used to 
probe the SZ effect, including the CSO (Caltech 
submm observatory) with BOLOCAM on Mauna 
Kea, ACBAR at the South Pole, the LMT (large 
mm telescope) in Mexico, and the LDB BLAST. 

A secondary anisotropy not associated with the 
CMB photons is stellar and accretion disk radi- 
ation reprocessed by dust into the infrared and 
redshifted into the submm, leading to a Wien tail 
distortion of the CMB blackbody. This has been 
found in the COBE FIRAS and shorter wave- 
length DIRBE data, with energy at a level about 
twice that in optical light. Anisotropies from dust 
emission from high redshift galaxies are being tar- 
geted by the JCMT with the SCUBA bolome- 
ter array, the CSO (soon with BOLOCAM), the 
OVRO mm interferometer, the SMA (submm ar- 
ray) on Mauna Kea, the LMT, the ambitious 
US/ESO ALMA array in Chile, the LDB BLAST, 
and ESA's FIRST satellite. About 50% of the 
submm background has so far been identified with 
sources that SCUBA has found. 

Ci's from nonlinearities in the gas at high red- 
shift are concentrated at high £, but for most vi- 
able models are expected to be a small contam- 
inant. Similarly, Thomson scattering from gas 
in moving clusters also has a small effect on Ci 
(although it should be measurable in individual 
clusters, e.g., with BOLOCAM and with Planck). 
The effect of lensing is to smooth slightly the 
acoustic peaks and troughs of Fig. and induce 
small scale non- Gaussian effects. 

1.3. Boomerang & Maxima 

Boomerang carried a 1.2m telescope with 16 
bolometers cooled to 0.3K in the focal plane aloft 
from McMurdo Bay in Antarctica in late Decem- 
ber 1998, circled the Pole for 10.6 days and landed 
just 50 km from the launch site, only slightly 
damaged. In @], maps at 90, 150 and 220 GHz 
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showed the same spatial features and the inten 
sities were shown to fall precisely on the CMI 
blackbody curve. The fourth frequency chan 
nel at 400 GHz is dust-dominated. Fig. |l| show 
the 150 GHz map derived using only one of th 
16 bolometers. Although Boomerang altogethe 
probed 1800 square degrees, only the region in th 
rectangle covering 440 square degrees was usee 
in the analysis described in |]]|| and this paper 
Fig. [j] also shows the 124 square degree regio: 
of the sky (in the Northern Hemisphere) tha 
Maxima-1 probed. Though Maxima was not ai 
LDB, it did so well because its bolometers were 
cooled even more than Boomerang's, leading to 
higher sensitivity per unit observing time; fur- 
ther, all frequency channels were used in creating 
its map. 

Analyzing Boomerang and other experiments 
involves a pipeline that takes (1) the timestream 
in each of the bolometer channels coming from 
the balloon plus information on where it is point- 
ing and turns it into (2) spatial maps for each 
frequency characterized by average temperature 
fluctuation values in each pixel (Fig. |l|) and a 
pixel-pixel correlation matrix characterizing the 
noise, from which various statistical quantities are 
derived, in particular (3) the temperature power 
spectrum as a function of multipole (Fig. ||), 
grouped into bands, and two band-band error ma- 
trices which together determine the full likelihood 
distribution of the bandpowers pc[ . Fundamen- 
tal to the first step is the extraction of the sky 
signal from the noise, using the only information 
we have, the pointing matrix mapping a bit in 
time onto a pixel position on the sky. 

There is generally another step in between (2) 
and (3), namely separating the multifrequency 
spatial maps into the physical components on the 
sky: the primary CMB, the thermal and kine- 
matic Sunyaev-Zeldovich effects, the dust, syn- 
chrotron and bremsstrahlung Galactic signals, 
the extragalactic radio and submillimetre sources. 
The strong agreement among the Boomerang 
maps indicates that to first order we can ignore 
this step, but it has to be taken into account as 
the precision increases. The Fig. || map is con- 
sistent with a Gaussian distribution, thus fully 
characterized by just the power spectrum. Higher 
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Figure 3. 2-sigma likelihood contours for the 
dark matter density and dark energy density (as- 
suming wq — — 1, i.e., a cosmological constant) 
are shown when progressively more CMB data are 
added to the LSS data. The smallest 2-sigma re- 
gion shows what happens when SNIa data is added to 
Boomerang+Maxima+prior-CMB+LSS. In all cases, 
the weak-H+age prior probability is adopted. Other 
variables have been marginalized, as in Table [lj 



order (concentration) statistics (3,4-point func- 
tions, etc.) tell us of non-Gaussian aspects, nec- 
essarily expected from the Galactic foreground 
and extragalactic source signals, but possible even 
in the early Universe fluctuations. For exam- 
ple, though non-Gaussianity occurs only in the 
more baroque inflation models of quantum noise, 
it is a necessary outcome of defect-driven models 
of structure formation. (Peaks compatible with 
Fig. ^ do not appear in non-baroque defect mod- 
els, which now appear unlikely.) 

1.4. Parameters of Structure Formation 

For this paper, we adopt the restricted set of 
7 cosmological parameters used in Among 
these, we have 2 characterizing the early uni- 
verse primordial power spectrum of gravitational 
potential fluctuations $, one giving the over- 
all power spectrum amplitude V^(k n ), and one 
defining the shape, a spectral tilt n s (k n ) = 
1 + <iln'P<j,/(ilnfc, at some (comoving) normaliza- 
tion wavenumber k n . We really need another 2, 
T^cwikn) and rit(k n ), associated with the gravi- 
tational wave component. In inflation, the ampli- 
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tude ratio is related to nt to lowest order, with 
0(n s — n t ) corrections at higher order, e.g., ||. 
There are also useful limiting cases for the n s — n t 
relation. However, as one allows the baroqueness 
of the inflation models to increase, one can en- 
tertain essentially any power spectrum (fully k- 
dependent n s {k) and nt(fc)) if one is artful enough 
in designing inflaton potential surfaces. As well, 
one can have more types of modes present, e.g., 
scalar isocurvature modes (Vi S {k n ),ni S (k)) in ad- 
dition to, or in place of, the scalar curvature 
modes (Vq,(k n ),n s (k)). However, our philoso- 
phy is consider minimal models first, then see 
how progressive relaxation of the constraints on 
the inflation models, at the expense of increas- 
ing baroqueness, causes the parameter errors to 
open up. For example, with COBE-DMR and 
Boomerang, we can probe the GW contribution, 
but the data are not powerful enough to deter- 
mine much. Planck can in principle probe the 
gravity wave contribution reasonably well. 

We use another 2 parameters to characterize 
the transport of the radiation through the era 
of photon decoupling, which is sensitive to the 
physical density of the various species of parti- 
cles present then, ojj = f2jh 2 . We really need 4: 
tOb for the baryons, Lu cc i m for the cold dark mat- 
ter, cohdm for the hot dark matter (massive but 
light neutrinos), and uj er for the relativistic par- 
ticles present at that time (photons, very light 
neutrinos, and possibly weakly interacting prod- 
ucts of late time particle decays) . For simplicity, 
though, we restrict ourselves to the conventional 3 
species of relativistic neutrinos plus photons, with 
u> er therefore fixed by the CMB temperature and 
the relationship between the neutrino and photon 
temperatures determined by the extra photon en- 
tropy accompanying e + e~ annihilation. Rather 
than be exhaustive when the data does not war- 
rant it, the effect of massive neutrinos is shown 
in Fig. ^ for a few values of ujhdm/^m, where 

l^>m = Uhdm + ^cdm + &b- 

The angular diameter distance relation adds 
dependence upon (1 — Q to t) and Q\ as well as 
on fi TO = u} m /h 2 |l2[ . This defines a functional 
relationship between Qtot and f^A, a degeneracy 
|TT| that would be exact except for the integrated 
Sachs- Wolfe effect, associated with the change of 
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Figure 4. 2-sigma likelihood contours in the (uj c dm + 
uJhdm) — plane are shown as the fraction of mat- 
ter in massive neutrinos is increased from up to 0.3, 
for Boomerang+Maxima+prior-CMB+LSS when the 
weak-H+age + Qtot = 1 prior probability is adopted. 
The age restriction used here, 11 < age < 16 Gyr, is 
slightly stronger than the > 10 Gyr weak-age prior 
adopted for the Table | limits. The CMB data 
alone cannot now discriminate among the cases - 
the left plot contrasts the large Ci variation as n s 
is varied over the tiny one as Qhdm is varied from 
to 0.95. When LSS is added, the contours shift in 
the expected direction, with more (cold+hot) dark 
matter preferred, at the expense of less dark en- 
ergy. Nonetheless, significant dark energy is appar- 
ently required. Note that the simplest interpreta- 
tion of the superKamiokande data on atmospheric 
is that f2„ T ~ 0.001, about the energy density 
of stars in the universe, which implies a cosmologi- 
cally negligible effect. Degeneracy between e.g., 
and v r could lead to the cosmologically very inter- 
esting Q. m v = Qhdm ~ 0.1, although the coincidence 
of closely related energy densities for baryons, CDM, 
HDM and dark energy required would be amazing. 
The uJb and n s likelihood curves are essentially inde- 
pendent of r2 m „/fi m for fixed f2 m . 
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$ with time if f^A or Qk is nonzero. (If <J> van- 
ishes, the energy of photons coming into potential 
wells is the same as that coming out, and there is 
no net impact of the rippled light cone upon the 
observed AT.) 

Our 7th parameter is an astrophysical one, the 
Compton "optical depth" tq from a reionization 
redshift z re h to the present. It lowers Cg by 
exp(— 2tc) at the high ^'s probed by Boomerang. 
For typical models of hierarchical structure for- 
mation, we expect tq < 0.2. It is partly degen- 
erate with tr 8 and cannot be determined at this 
precision by CMB data now. 

The LSS also depends upon our parameter set: 
the most important combination is the wavenum- 
ber of the horizon when the energy density in 
relativistic particles equals the energy density in 
nonrelativistic particles: k~^ eq w 5r _1 h~ 1 Mpc, 

where T w O m hf2 el V 2 . Instead of V$>{k n ) for the 
amplitude parameter, we often use Cio at £ = 10 
for CMB only, and a\ when LSS is added; of is a 
bandpower for density fluctuations at 8h _1 Mpc, 
a scale associated with rare clusters of galaxies. 
When LSS is considered in this paper, it refers to 
constraints on T + (n s — l)/2 and In of that are 
obtained by comparison with the data on galaxy 
clustering and cluster abundances 0. 

When we allow for freedom in Lu er , the 
abundance of primordial helium, tilts of tilts 
{d n {s,is,t}{k n ) I d\n.k, ...) for 3 types of perturba- 
tions, the parameter count would be 17, and 
many more if we open up full theoretical freedom 
in spectral shapes. However, as we shall see, as of 
now only 3 or 4 combinations can be determined 
with 10% accuracy with the CMB. 

Table [I] shows there are strong detections with 
only the CMB data for f2 tot , coi, and n s in the min- 
imal inflation-based 7 parameter set. The ranges 
quoted are Bayesian 50% values and the errors 
are 1-sigma, obtained after projecting (marginal- 
izing) over all other parameters. With Maxima, 
oJcdm begins to localize, but much more so when 
LSS information is added. Indeed, even with just 
the COBE-DMR+LSS data, the 2-sigma contours 
shown in Fig. || are already localized in oj c dm- 
That Q\ is not well determined is a manifestation 
of the fltot-^A near-degeneracy discussed above, 



which is broken when LSS is added because the 
CMB-normalized as is quite different for open cf. 
A models. Supernova at high redshift give com- 
plementary information to the CMB, but with 
CMB+LSS (and the inflation-based paradigm) 
we do not need it: the CMB+SN1 and CMB+LSS 
numbers are quite compatible. In our space, the 
Hubble parameter, h = (X^^jh 2 )) 1 / 2 , and the 
age of the Universe, to, are derived functions of 
the fljh 2 : representative values are given in the 
Table caption. 

Since this is a conference devoted to neutrinos, 
Fig. |4| is included to show what happens as we let 
the fraction of the matter in massive neutrinos 
vary, from to 0.3. Until Planck precision, the 
CMB data by itself will not be able to strongly 
discriminate this ratio. Adding HDM does have a 
strong impact on the CMB-normalized as and the 
shape of the density power spectrum (effective T 
parameter), both of which mean that adding some 
HDM to CDM is strongly preferred in the absence 
of f^A- However a nonzero (though smaller) Q\ 
is still preferred. 

We can also forecast dramatically improved 
precision with further analysis of Boomerang and 
Maxima, other future LDBs, MAP and Planck. 
Because there are correlations among the phys- 
ical variables we wish to determine, including 
a number of near-degeneracies beyond that for 
Qtot-^A [0 ! it is useful to disentangle them, by 
making combinations which diagonalize the er- 
ror correlation matrix, "parameter eigenmodes" 
]qjll|]. For this exercise, we will add uihdm and 
n t to our parameter mix, making 9. (The ra- 
tio VGw{k n ) /V^ikn) is treated as fixed by n f , a 
reasonably accurate inflation theory result.) The 
forecast for Boomerang based on the 440 sq. deg. 
patch with a single 150 GHz bolometer used in the 
published data is 3 out of 9 linear combinations 
should be determined to ±0.1 accuracy. This is 
indeed what we get in the full analysis of the Ta- 
ble and Boomerang+DMR only. If 4 of the 6 
150 GHz channels are used and the region is dou- 
bled in size, we predict 4/9 could be determined 
to ±0.1 accuracy. And if the optimistic case for 
all the proposed LDBs is assumed, 6/9 param- 
eter combinations could be determined to ±0.1 
accuracy, 2/9 to ±0.01 accuracy. The situation 
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Table 1 

Cosmological parameter values and their 1-sigma errors are shown, determined after marginalizing over the other 
6 cosmological and 4 + experimental parameters, for B98+Maxima-I+prior-CMB and the weak prior used in 
(0.45 < h < 0.9, age > 10 Gyr). The LSS prior was also designed to be weak. The detections are clearly 
very stable if extra "prior" probabilities for LSS and SN1 are included. (Indeed, they are stable to inclusion of 
stronger priors -except if the BBN-derived 0.019 ± 0.002 is imposed §.) Similar tables for B98+DMR are given 
in Q and for B98+MAXIMA-I+DMR in If Q.tot is varied, parameters derived from our basic 7 come out to 
be: age=13.2 ± 1.3 Gyr, h = 0.70 ± 0.09, D. m = 0.35 ± .06, Q b = 0.065 ± .02. Restriction to Q, tot = 1 yields: 
age=11.6 ± 0.4 Gyr, h = 0.80 ± .04, Q m = 0.31 ± .03, Q b = 0.046 ± .005. 
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improves for the satellite experiments: for MAP, 
we forecast 6/9 combos to ±0.1 accuracy, 3/9 to 
±0.01 accuracy; for Planck, 7/9 to ±0.1 accuracy, 
5/9 to ±0.01 accuracy. While we can expect sys- 
tematic errors to loom as the real arbiter of accu- 
racy, the clear forecast is for a very rosy decade of 
high precision CMB cosmology that we are now 
fully into. 
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